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We investigate the behavior of a harmonically trapped system consisting of an impurity in a dilute 
ideal Bose gas after the boson-impurity interaction is suddenly switched on. As theoretical frame¬ 
work, we use a field theory approach in the space-time domain within the T-matrix approximation. 

We establish the form of the corresponding T-matrix and address the dynamical properties of the 
system. As a numerical application, we consider a simple system of a weakly interacting impurity 
in one dimension where the interaction leads to oscillations of the impurity density. Moreover, we 
show that the amplitude of the oscillations can be driven by periodically switching the interaction 
on and off. 


PACS numbers: 67.85.Pq. 

I. INTRODUCTION 

Understanding the behavior of an impurity in a host 
medium is a fundamental question in quantum many- 
body physics. On the one hand, this is motivated by 
the possibility to use an impurity particle as a probe in 
experiments. On the other hand, such systems help to 
gain insights and build analytical and numerical tools 
to approach more complicated impure systems theoret¬ 
ically. A standard example is an electron placed in an 
ionic crystal, the solution to which can be approached by 
introducing the concept of the polaron - a quasiparticle 
that can be visualized as an electron dressed by collec¬ 
tive excitations in the medium. This concept was first 
proposed in 1933 by Lev Landau [1] and has been widely 
used ever since [2-4]. One particularly interesting ques¬ 
tion in the study of impure systems is the time dynamics 
of an impurity in a medium. This question was mainly 
addressed for charged impurities in polar lattices or or¬ 
ganic crystals [5-13] and still represents a challenge for 
the theoretical description. 

Modern experiments with cold atomic gases produce 
systems where the density of immersed atoms is much 
smaller than the density of the host medium which can 
consist of fermionic [14, 15] or bosonic particles [16, 17]. 
The scattering length between host and impurity par¬ 
ticles can be tuned using Feshbach resonances and the 
spatial geometry of the system can be designed using ex¬ 
ternally modulated traps [18-22]. It is also worth to note 
hybrid set-ups which contain a single ion in a neutral sys¬ 
tem [23, 24] or an electron in a Bose-Einstein condensate 
(BEG) [25], which leads to a different shape of the inter¬ 
action potential in the system. 

The experimental progress in the realization of impure 
cold atomic systems allows to investigate the associated 
quasiparticles that appear due to the interaction between 
an impurity and the medium. In the literature these 
quasiparticles in a Eermionic (Bosonic) medium are of¬ 
ten referred to as Eermi (Bose) polarons in analogy with 
Landau’s concept. Experimental realization of degener¬ 
ate imbalanced Eermi gases and subsequent observation 


of a polaron in a Eermi sea [14], drove a wave of new 
theoretical activity in the Eermi polaron problem where 
a number of reliable results has been produced using dif¬ 
ferent Monte Carlo algorithms [26-28], variational wave 
functions [29, 30] and T-matrix approaches [31]. Eor a 
review of these activities, see Ref. [32] and references 
therein. 

The experimental and theoretical study of an neutral 
impurity in a system of cold bosons has not yet been pur¬ 
sued as extensively as in the fermionic case. The theoret¬ 
ical approaches were motivated by systems of superfluid 
^He doped with ^He, neutrons, etc. and were formulated 
for homogeneous systems [33, 34]. They have usually 
been applied to weakly-interacting dilute systems. Low- 
density cold atomic gases offer a new playground for the¬ 
oretical studies. In these systems, the impurity-boson in¬ 
teraction can often be regarded of zero-range with a tun¬ 
able coupling constant. The temperature of the medium 
can be very low and as a first approximation it can be 
put to zero. 

If the interaction in the system is weak, the usual ap¬ 
proach to calculate observables relies on mean-field equa¬ 
tions, e.g. the Gross-Pitaevski equation. These equa¬ 
tions have allowed the study of the motion of an impurity 
through a BEG in different spatial dimensions [35] and 
self-localization of an impurity in a trapped BEG [36]. 
Although the mean-field equations are limited in interac¬ 
tion strength they usually offer a natural way to include 
the effects of the trap which is present in experiments. 
To study regimes where the boson-impurity interaction is 
strong, one can use variational wave functions to extract 
an upper bound for the energy as well as the structure of 
the system [37-41]. 

In cases where the Bogoliubov approximation for the 
condensate mode is applicable, a commonly used ap¬ 
proach is to map the original Hamiltonian onto the 
Erohlich Hamiltonian [2]. Eor cold atomic gases this 
Hamiltonian was approached using the Eeynman path in¬ 
tegral formulation [42], diagrammatic Monte Garlo [43], 
the renormalization group approach [44], and using a 
variational wave functions in the form of correlated Gaus- 
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sians [45]. It was also shown that using the T-matrix 
approximation, the problem can be approached almost 
analytically to some extent [46]. For time-independent 
homogeneous systems, these methods are able to predict 
the effective mass of the polaron for weak and strong 
coupling regimes. 

The present work continues the study of an neutral 
impurity immersed in a Bose gas at zero temperature. 
Using the field-theoretical approach and the T-matrix ap¬ 
proximation, we investigate the real-time dynamics of an 
impurity in a trapped ideal Bose gas after the impurity- 
boson interaction is suddenly turned on. We first work 
out the form of the in-medium T-matrix. This knowl¬ 
edge allows us to investigate the dynamical properties 
of the system. As a numerical application, we consider 
a weakly-interacting one-dimensional system. After the 
interaction is turned on, the impurity oscillates in the 
trap. We also show that by turning the interaction on 
and off periodically one can drive the amplitude of the 
oscillations. 

The paper is organized as follows: in Section II we 
review the basics [3, 47] of the field theory approach 
to many-body problems in the space-time domain with 
a particular emphasis on harmonically trapped systems 
and T-matrix approximation. In Sections III and IV, we 
derive the T-matrix for systems of different spatial di¬ 
mension analytically and discuss our approach to dynam¬ 
ical properties. This is illustrated by a numerical appli¬ 
cation to the non-trivial example of a weakly-interacting 
one-dimensional system. Our conclusions and a discus¬ 
sion of possible applications and extensions of this work 
are given in Section V. Finally, we include four appen¬ 
dices with technical details and explicit derivations used 
in the main text. 


II. FORMULATION OF THE PROBLEM 


Hamiltonian. We consider a system with one impurity 
particle of mass mj and N non-interacting bosons of mass 
in I) spatial dimensions. The Hamiltonian for this 
system consists of three pieces 


H = ( 1 ) 


where Hb describes the system of bosons, Hj is the 
Hamiltonian for the impurity, and W is the impurity- 
boson interaction. The ideal, harmonically trapped Bose 
gas is described by the Hamiltonian 


Hb 



2mB dx‘^ 


2 


0(x)dx, 


( 2 ) 


whereas the Hamiltonian for the impurity has the form 


Hi 



h? miQ?x‘^ 

2m/ dx‘^ ^ 2 


'0(x)dx, 


( 3 ) 


where x is a coordinate (scalar in one spatial dimension or 
vector in two or three spatial dimensions) and are 


the field operators that annihilate a boson (impurity). 
These operators are defined in a standard way, for in¬ 
stance (j){x) = where bi is the annihilation 

operator in seconcTquantization and (x) is the normal¬ 
ized one-body eigenfunction of Hb- For clarity, the field 
operators will be denoted with lower-case letters while 
the wave functions will be written with corresponding 
capital letters. 

The interaction between the impurity and the bosons 
is switched on at time t = 0. It can be written as 

W = 0{t) J (1)^{xi)'lp^{X 2 )V{xi - X 2 )(l){xi)'ip{x 2 )dxidx 2 , 

( 4 ) 

where V is the impurity-boson interaction potential and 
0{t) is the Heaviside step function.^ Initially the system 
is in the non-interacting ground state, but the dynam¬ 
ics becomes non-trivial after the interaction is turned 
on at t = 0. Note, that introducing a trap makes the 
calculations much more inovolved compared to the ho¬ 
mogeneous case. However, some aspects of the problem 
become more straightforward. For example to describe 
the homogeneous case one often works in a finite space 
parametrized by some length, L, and then takes the limit 
N{L) oo, NjL^ const. In the trapped case the sys¬ 
tem is localised and taking this limit is not required. 

Wave function for impurity. To describe dynamics of 
the impurity after for t > 0, we focus on the following 
object 



N 

F{y\Zi,...,ZN\t) W*o{zi)dzi, 

i=l 


( 5 ) 


where F{y; zi,..., znH) is the wave function of the system 
of one impurity with coordinate y and N bosons with co¬ 
ordinates 2 : at time t. In the present paper, we consider 
a very dilute Bose gas for which the ladder approxima¬ 
tion discussed below holds. In this approximation the 
density of the bosons is not changed drastically, so the 
object I{yH) can be considered as the wave function of 
the impurity at time t, even though the impurity alone 
cannot be treated as an isolated system. If we know the 
wave function for the impurity at t = 0, then I{yH) foi* 
t > 0 can be obtained as 


j dxGi{x,0,y,t)I{x,0), 


where G/ is the Green’s function for the impurity 


( 6 ) 


W{oo)) 


^ Note that the effective single-channel description implied by the 
interaction in Eq. (4) can be used for broad Feshbach resonances 
[18]. 
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defined using field operators in the interaction picture, 
i.e. 2 pi{x,t) = . Here T is the stan¬ 

dard time ordering operator (larger times are always to 

the left); and S{t) = T exp ^ Wi(t')dt'^. The an¬ 
gle brackets in Eq. (7) indicate averaging over the ideal 
Bose gas. As we consider only one impurity in a non¬ 
interacting Bose gas, (5'(oo)) = 1. In the diagram tech¬ 
nique this is very easily seen, since with one impurity and 
no boson-boson interactions only connected diagrams are 
possible. It is obvious that G/(x, t,x', t') should depend 
only on t' — t, since the presented formulation can also 
correspond to a system where the boson-impurity inter¬ 
action is constant in time and the impurity is immersed 
in the medium at time t. 


Bogoliubov approximation. The non-interacting Bose 
gas at zero temperature has all particles in the ground 
state of the one-body Hamiltonian. To deal with macro¬ 
scopic number of particles occupying one level, we uti¬ 
lize the so-called Bogoliubov prescription. First we 
write the field operators for the particles in the conden¬ 
sate and particles depleted from the condensate explic¬ 
itly, i.e. (j){x,t) = 6o^o(^) d-(l)'i{x,t), where = 

^iHstlh write this expres¬ 


sion down we assume that the energy of the condensate is 
zero, which for a harmonic trap technically corresponds 
to subtracting a constant term DM112 from Hb and Hj. 
After the operators for the condensate and the depleted 
bosons are separated, we obtain four interaction terms in 
Eq. (4). In these four terms, assuming that the number 
of depleted bosons is small, we perform the substitution 
6o — 6o — ^/No^ where Nq is the number of particles 
in the condensate. This number should be treated as a 
variable that minimizes the total energy. We want to 
investigate the impurity in a dilute gas at zero tempera¬ 
ture assuming that the density of the condensate is not 
appreciably changed due to the impurity. This is done in 
the so-called ladder approximation with only one particle 
depleted from the condensate at any moment of time. So 
if the the total number of bosons is large, i.e. A" ^ I, we 


may substitute the operators bo and 6 q with \^N. 

Ladder approximation. We consider a dilute gas, such 
that I, where R defines the relevant interac¬ 

tion range. The density n is defined as the peak density of 
the BEC^ in the trap. If {n^^^R) I, the system can be 
approached using the ladder approximation where only 
one boson is depleted at any moment of time 

Gi{x,t,x',t') = + N ( Axidyidti^Q{yi) 

i=l,2 


^Q{y2)G°{x,t,Xi,ti)r{xi,yi,ti,X2,y2,t2)Gi{X2,t2,x',t'), 



FIG. I. Diagrams representing main equations of the lad¬ 
der approximation discussed in the text, (a) Eq. (8) for the 
impurity Green’s function, (b) Eq. (9) for the in-medium T- 
matrix. Solid thick (thin) lines stand for the impurity Green’s 
function, Gi (G/). The filled boxes represent the in-medium 
T-matrix, F. The wavy lines denote the interaction potential 
V. The zig-zag lines depicts the condensate particles, and the 
dashed line is the free Green’s function for a boson depleted 
from the condensate, G%. 


where F is the in-medium T-matrix discussed below 
and G^ is the free impurity propagator determined from 
Eq. (7) with W = 0. Note that in this paper the inte¬ 
gration limits Too are implied if nothing else is specified. 
It is worth to note that the integration over time vari¬ 
ables is effectively always done over a finite interval, since 
both G 5 (xi,ti,X 2 ,t 2 ) and r(xi,^i,G,X 2 ,^ 2 ,^ 2 ) vanish 
for G < ti. This observation follows since the impurity 
cannot propagate backwards in time. The diagrammatic 
representation of Eq. (8) is presented in Fig. I. The in¬ 
medium T-matrix, F, is written as an integral equation. 


F(xi,^i,ti,X 2 ,^ 2 ,G) =^(^1 - X 2 )S{yi - y 2 )S{ti - G)E(xi - yi) 

+^V{xi-yi) j G°{xi,ti,y,T)G%{yi,ti,v',TX{y,v',T,X2,y2,t2)<^y<^y'^T , (9) 


^ Note that, for simplicity, we generically use the term BEC for a 
system with microscopic occupation number of the ground state. 
This system, however, has different properties in different spatial 
dimensions [48]. 


where G^ is the free Green’s function for depleted bosons 
(cf. Fig. I). For the homogeneous case, the ladder ap¬ 
proximation in the momentum-frequency domain was 
shown to capture the qualitative behavior of the impurity 
in the Bose-polaron problem [46] . 

Free Greenes funetions in a harmonie trap. In a ID 
harmonic trap, the free Green’s function for t > 0 can be 
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written as 


G°j{x,0,x',t) = , (10) 

n=0 


with Cn = hfln (recall that we subtract the zero point 
motion). The corresponding one-body wave function can 
be always chosen as a real function, i.e. 




.1/4 

e (miil \ ' ^ 



( 11 ) 


where Hn is a Hermite Polynomial of order n. From 
Eq. (10), we can derive a useful identity that connects 
different time scales: 


G°{x,0,x',t') =i j G°{x,0,y,t)G°{y,t,x\t')dy. (12) 

Since a D-dimensional harmonic trap decouples in Carte¬ 
sian coordinates, the corresponding Green’s function is 
just a product of ID Green’s functions from Eq. (10). 
This product can be rewritten (e.g. using the Mehler for¬ 
mula) in the following form 


G5(x, 0, x', t) = —ie* 2 


. D^t f miCl 


2i7rhsm{flt) 


exp 


imiQ 
2hsm{Qt) 


(^{x^ -h x'^) cos(flt) — 2xx'^ 


(13) 


where the terms of the form xx' should be understood 
as dot products. The ID Green’s function for the bosons 
depleted from the condensate is written as 


oo 

GB{x,0,x',t) = -* E , (14) 

n=l 

where 4>^(x) is given by Eq. (11) with rriB instead of 
mi. We see that can be written in a form similar 
to Eq. (13) up to the extra term i4>o(x)4>o(x') since in 
Eq. (14) the sum starts from one. This extra term is 
constant in a homogeneous case and can be neglected. It 
can also be disregarded if one considers only the leading 
order perturbation from the weak interacting potentials, 
where one neglects higher order corrections in T. More¬ 
over, this term also gives higher order contributions in the 
full Green’s functions if one considers short time scales, 
QT 1, since in this case the leading effect arises from 
the propagation in free space. To get rid of this term 
completely one can redefine G% to also include n = 0 
term. Physically, this corresponds to depleted particles 


remaining in the ground state of the trap which prop¬ 
agate with a phase different from the condensate. In 
our discussion below, we adopt this reformulation of the 
Green’s function, which now reads 


G^(x,0,x',t) = —ie^ 


pnt f ms^ 


27rhi sin(flt) 


X D/2 

) " 


exp 


imB^ 


(^{x‘^ -h x'^) cos(flt) — 2xx'^ 


2hsm{flt) j yj (5) 

These expressions can be easily written in the 
momentum-time domain. Eor simplicity in the deriva¬ 
tions below we will use wave vectors instead of momenta. 
We define the corresponding Green’s function as 


G(pi,ti,p2,t2) = [ G(xi,ti,X2,t2)e*^^^^ ^^^P^dxidx2. 

(16) 

The signs in the argument of the exponent are deter¬ 
mined from the Eourier transform of the creation (anni¬ 
hilation) operators in the definition of the Green’s func¬ 
tion. The Green’s function in momentum space for the 
impurity is 


G°j{k,0,k',t) = (_ 

\ir 


27rh 


D/2 


exp 


im/fl sin(flt) ^ 

- ' {{k^-\-k'^)co^{VLt) — 2kk'\ , (17) 

_2imiVL^m{VLt) ^ ^ ^ ^ G 


and for the bosonic particles the mass m b should be used 
instead of mi. 

Interaction. We would like to investigate universal be¬ 
havior of the system so any particular form of the poten¬ 
tial will work. In this paper, we employ a simple Gaus¬ 
sian potential. 


V{x) = 9(^-^'' (18) 

where g and A are parameters that determine the inter¬ 
action. Note that for A —^ oc the interaction becomes of 
zero range. Eor the derivations, we will use this limit as 
we cannot describe short-distance physics anyway. In two 
and three spatial dimensions one has to relate the cou¬ 
pling strength g to the scattering length before taking the 
zero-range limit in order to describe physical observables 
(see Appendix A). 


III. T-MATRIX 

Let us see what happens if the limit A ^ oo is taken 
right away which yields V{x) = gS{x). We know that in 
free space T corresponds to a dimer propagator, so we 
look for T in the following form 
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' 5{x, - y,)5{x, - f ^ 


2'Khi sm{Qt) J 


( 19 ) 


where t > 0. The function /(t) describes ’formation’ and 
’decay’ of this virtual (coherent) dimer. The form of T in 
Eq. (19) is obvious by noticing that the frequency is the 
same for both particles. Therefore the product 
in Eq. (9) can be interpreted as a dimer propagator and 
Eq. (12) can be applied. Eor more complicated trapping 
potentials the form of T is less obvious. 

It is worth to note that f{t) determines the time scale 
at which the diagrams with the T-matrix become im¬ 
portant when compared to the free propagation. Eor in¬ 
stance in the example we consider in the next section, we 
use f{t) = gS{t), which means that for t h/{gn) the 
particle can be well described using only the free propa¬ 
gator. In the same way by comparing the diagrams that 
are not in the ladder approximation with the diagramms 
that are included, the function f{t) can be used to vali¬ 
date the use of the ladder approximation in a dilute gas. 

The equation for / can be obtained from Eq. (9) 


fit) = gKt) - j 



where fi = miVfiB j{jni + m^) is the reduced mass and 
the lower integration limit A ^ 1 /denotes that the in¬ 
tegral should be regularized as otherwise it appears to be 
divergent \i D > 1. The need to introduce this regulator 
is related to the use of zero-range interactions for which 
^ 0 if the scattering length is fixed (see Appendix A), 

so the limit of A ^ oc should be taken properly. This 
procedure yields finite quantities expressed through the 
scattering length as shown in Appendix B for 2D and 
3D. Equation (20) is already suitable for numerical cal¬ 
culations. However, since the equation is linear in / we 
perform a Eourier transform to obtain further analytic 
insight: 


fit) = ^ j (21) 

which leads to the equation for /(cj). 


/(w) = g- 


wfii^) 


—i 

h 

\2TrhJ Jx 



( 22 ) 


where we used that f{t) is zero for t < 0. To make 
sense of the integral in the expression, we imply that uj 
has a vanishing positive imaginary part, uj cj -h iO. 
It is worth noting that one can restrict the integration 
limits in eq. (22)using periodicity of the sin(flt) function 
which might be useful, for example, for regularization of 


the integral. Now we examine this equation in different 
spatial dimensions. 

One-dimensional ease. We first consider the ID case 
where we have 


/(w) = g- 


igfi^^) 

h 



dr. (23) 


Eor cold atomic gases experiments in highly elongated 
’’cigar”-shaped optical traps it was shown [49] that 


2h?a I 


(24) 


where a is the 3D scattering length, L is the harmonic 
oscillator length in the squeezed direction, and G = 
|C(I/2)| = 1.46..., where C is the Riemann zeta function. 
The integral in Eq. (23) can be easily taken such that 


/(‘^) 


g + 


gfi^) / M p 
h V 


/I CJ 
V2 ’m 



where B{x^y) is the standard beta function. Eor the 
sake of argument here and below we will assume that 
^ 0, which corresponds to considering short time 
scales fit 1. It is worth to note that the results de¬ 
termined in this way should have the same form as for 
the homogeneous case, where = 0, so fit ^ I is triv¬ 
ially satisfied. This assumption allows us to demonstrate 
the main ideas in a relatively simple manner. Eor practi¬ 
cal purposes /(t) can be easily found numerically either 
by using the inverse Eourier transform or by solving the 
corresponding equation directly in the time domain. 

Assuming short times leads to 


I 

( 26 ) 

where the superscript S means that fit ^ I. The in¬ 
verse transformation can be easily performed for large 
and small values of g. Indeed if ^ 0 we obtain 
y (w) ~ g(l - f\/^^), which corresponds to f^it) = 

gS(t) — These terms as well as higher or- 

der contributions can be derived directly in the time do¬ 
main (see Appendix C). Eor repulsive interaction, ^ > 0, 
one can easily perform inverse Eourier transform 





\/2/j7r 


+ ^erfcx 


V ) \ ’ 
(27) 
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where erfcx(x) is the scaled complementary error func¬ 
tion. This form corresponds to the limit of the infinite 
sum derived in Appendix C. For strong repulsive interac¬ 
tion, i.e. ^ ^ oo, for very short times, i.e. t 

function f^{t) is the same as for the weakly-interacting 
case discussed above. For larger times we derive 


fit) ^ h 


i — 1 



(28) 


Note that, this behavior of f{t) leads to non-analytic 
behavior of Gj for infinite repulsion. We understand that 
this shows strong correlations taking place in ID. These 
correlations lead to depletion of the condensate which is 
not properly taken into account in the present approach. 

Three-dimensional ease. After proper regularization 
(see Appendix B) we derive in 3D 


fif 


1 


2nh?a 



VcJ 

\/2nh 


(29) 


Notice that this equation with uj ^ uu— 2h{rni+mB) 
compared with the self energy derived in Ref. [46] for the 
homogeneous case, i.e. with ^4 = 0 , which corresponds 
to the pair propagator in free space. It is also interesting 
to note that the functional form of Eq. (29) is similar to 
Eq. (26). Apparently the parametric dependence on g 
in ID is the same as dependence on — 1/a in 3D, which 
means, for instance, that the weakly interacting ID case 
has the same time dependence as the strongly interacting 
3D case. It also means that Eq. (27) with some minor 
changes in constants can be applied to find the inverse 
Eourier transform for a < 0 in 3D. Let us discuss the case 
of a < 0 now in more details. Apparently, if a —> 0 , f^{t) 
to leading order is a delta function. Obviously, this form 
is due to neglecting higher order diagrams such that the 
virtual dimer does not propagate in time. At unitarity 
when 1/a = 0 , the expression for f^{t) can also be found 
easily. We note that has a branch cut for nega¬ 
tive values of cj, where we take ^ 1/\Juj -h id as 

discussed above. This leads to vanishing /^(t) for nega¬ 
tive values of t. This expression can be easily transformed 
back to the time domain, where the following form should 
be used 


fit) = -il + i)VTT 


Vi’ 


(30) 


which up to the delta function mimics time dependence 
of f^{t) for weak interaction in ID. The case of a > 0 
is more involved, since in 3D it implies existence of both 
a bound molecular state and repulsive polaron branches, 
see e.g. Ref. [46]. We briefly address this problem for 
small positive a in the time domain in Appendix D. 

Two-dimensional ease. Einally, we briefly mention the 
2 D case. After regularization (see Appendix B), we ob¬ 
tain the following expression 


fiV 


2h?TT 1 

g iir — \r].{a2jjjaue^/h)^ 


(31) 


where 020 is the two-dimensional scattering length, and 
7 = 0.5772... is Euler’s constant. Again, we note that 

2 

this equation with uj ^ uj — 2h{rYii+mB) Vised to 

determine the self-energy without the trap. 


IV. DYNAMICS 

After the form of the T-matrix is established one can 
approach the wave function for the impurity, 

00 n 

0 = E n(^^?r)*G?/(t = 0), (32) 

n=0 i=0 

where operator notation is assumed. This equation is ob¬ 
tained by combining eq. ( 6 ) with the Green’s function 
from eq. ( 8 ) written as an iterative sum. In this expres¬ 
sion the integrals over the space variables can easily be 
taken, since the functions are just Gaussians and the in¬ 
tegrals over time variables can be performed numerically. 
Notice that higher order terms contribute to longer time 
evolution and are usually suppressed if one is interested 
in short-time dynamics. 

We turn our attention to a one-dimensional weakly- 
interacting system where the interaction is switched on 
at time t = 0 and the impurity starts to move in the trap. 
This simple system allows us to illustrate time dynamics 
almost analytically. In this paper, we are mostly focusing 
on the short-scale time dynamics, however, a quench of 
the coupling constant in the Hamiltonian can also be used 
to improve our understanding of non-equilibrium dynam¬ 
ics and relaxation processes in closed quantum systems, 
see e.g. [50, 51] and references therein. 

The motion of the impurity is described with the func¬ 
tion I{y,t) in Eq. (32) where for the sake of argument 
we assume that the system is weakly-interacting and 
f(t) gd{t). To find we can use Eq. (32) di¬ 

rectly. Another way is to notice that analog to the free 
Green’s function for the weakly-interacting case, we have 

Giix,t,x',t') = i J dyGiix,t,y,T)Giiy,T,x',t'), (33) 

which allows us to construct the infinitesimal operator of 
time evolution 

Gj(x, t, x', t') = G5(a:, t, x', t') + 

gN J dydT^l{y)G°ix,t,y,T)G°iy,T,x',t'), (34) 

and then apply it repeatedly to the initial wave func¬ 
tion. The integration over spatial coordinates can easily 
be taken, since only Gaussian functions enter. We are 
then left only with the integration over the time variables 
which we perform numerically. It is worth to note that 
Eq. (33) cannot be applied to more strongly interacting 
systems, where correlations between different times are 
more complicated. 
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We first discuss the homogeneous case, i.e. 0 = 0. In 
this case the infinitesimal time evolution operator will be 
written as 

G/(x, = (1 — — t)gn/h)G^j{x,t,x' (35) 

where n is the density of BEG. Notice that this operator 
does not conserve the normalization of the wave function 
in the order g^{t 'due to the infinitesimal character of 
the operator. After repeated use, we obtain the following 
Green’s function 

G/(x, t, x', t') = ~G^j{x^ t, x\ G), (36) 

which describes the shift of the energy due to interaction. 
This form of the Green’s function implies that if the im¬ 
purity is in an eigenstate, it will remain in the same state 
with some additional phase oscillations. In particular, 
if /(x,0) is an eigenstate, then I{x,t) = 0). 

This approach can be regarded as a mean-field treatment 
of the problem, i.e. we assume that the bosons form a 
homogeneous external field, which is not affected by the 
impurity. If the condensate density is not constant then 
these phase oscillations will be position-dependent, lead¬ 
ing to a probability current. One would expect this fiux 
initially to be proportional to the gradient of the density 
of the condensate multiplied by gP{x^ 0). At later times 
additional fluxes due to the external trap will appear. 
Since in the equations only the combination g^Q{x) en¬ 
ters, the same behavior emerges also if the condensate is 
homogeneous but the coupling constant is a function of 
the coordinates. 

To describe this system, we need to find G^(x, t, x', t') 
which now has a more complicated form compared to 
the homogeneous case. However, using simple numeri¬ 
cal integration routines the wave function /(^, t) still can 
be obtained. To simplify the notation, we use the units 
where h = mi = Cl = 1. This corresponds to measuring 
time in units of l/f^, energy in units of hCl^ and distances 
in units of ^yh/{mjCl). For our numerical calculations, 
we choose mj = ms and gN = 0.3. The corresponding 
density of the impurity as a function of time is shown in 
in Fig. 2. 

We see that the initial density starts to oscillate with 
almost complete revival after some time. This is eas¬ 
ily understood noticing that a weak external potential 
gN^Q{x) acting on the impurity would yield the same 
diagrams as we use now. Apparently, this external po¬ 
tential changes the eigenstates of the Hamiltonian such 
that the dynamics in the system should be described with 
the following Green’s function 

oo 

Gi{x,0,x',t') = (37) 

n=0 

where Sn is a new eigenenergy and F{x) is the corre¬ 
sponding eigenfunction. Using first order perturbation 



FIG. 2. (Color online) Density of the impurity, |/(x,t)|^ as a 
function of time, t, and position coordinate, x. The interpar¬ 
ticle interaction is assumed to be constant in time, gN = 0.3. 
t is in units of 1/Q, x is in units of y/KJJrniCi) and \I(x^t)\^ 
is in units of (miCl)/h. 

theory [52] , we approximate these quantities as 

£n-n + gNann, Onk = j ^'))(a:)$o{a;)^'fe(a;)da;, 

Fn{x) nix) +gNj2 ^kix). (38) 

kjtn 

Assuming /(x,0) = To(x), the wave function at later 
times up to the phase factor reads 

lix, t) = ^'o(^) FgNY, _ 

(39) 

We checked that this wave function reproduces the den¬ 
sity from Fig. 2 very accurately. Noticing that aon in 
Eq. (39) is non-zero only for even n due to the parity 
conservation, it becomes obvious that after t = ir the 
density almost completely returns to its initial profile. 
This revival is due to the equidistant energy spectrum of 
the harmonic trap so the period of oscillations is deter¬ 
mined by the oscillator frequency. The coupling constant 
and the density of bosons determine the amplitude of the 
oscillations. 

From previous studies in solid state physics [8, 11], it is 
known that the bare charge carrier evolves into a polaron 
by emitting phonons. The time dynamics of such a pro¬ 
cess strongly depends on the parameters of the problem. 
At the same time Fig. 2 indicates that the impurity does 
not evolve into a polaron at the time scales given by the 
trap. This result arises since our model forbids the deple¬ 
tion of the condensate in the final state. In the ideal case, 
this approximation should be justified. Indeed, there the 
excitation processes are of one-body nature and thus the 
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depletion processes at the presented time scales are sup¬ 
pressed by 1/\/N compared to non-depleting processes. 
To arrive at this conclusion, it is crucial that the interac¬ 
tion between the bosons is absent o r very we ak. In this 
context Very weak’ means NaBB\^^/{miQ)/L‘^ ^ 1 in 
the quasi-one-dimensional regime. If this condition is 
satisfied then the system behaves as an ideal gas [53]. 

Even though one can engineer effectively non¬ 
interacting systems in a laboratory by tuning the boson- 
boson scattering length [18], the question of what 

happens if the system posesses a weak boson-boson in¬ 
teraction is important. In this case, the boson-boson 
interaction term. 


1 

2 



ldX 2 (p'^ {xi)(p\x 2 )Vbb{xi - X2)(j>{xi)4>{x2), 


(40) 


should be added to the Hamiltonian. If this term is in¬ 
cluded, the excitation processes are of collective nature, 
see e.g. Ref. [53, 54]. Then one should take the deple¬ 
tion of the condensate explicitly into account, which will 
lead to the formation of a polaron by emitting collective 
excitations. The details of the formation, however, will 
depend strongly on the strength of the boson-boson in¬ 
teraction and are left for the future studies, i.e., in the 
present paper we always put V^h = 0. For relevant stud¬ 
ies with Vhh 7 ^ 0 in ID see Ref. [55] and an investigation 
of a three-body ensemble. Ref. [56]. 

Figure 2 graphically demonstrates that the system has 
some intrinsic time scale, so if one designs the coupling 
strength such that it depends on time it can give addi¬ 
tional tools to study dynamics. For instance if the inter¬ 
action is turned off after time t = tt the oscillations of 
the density would be reduced, whereas if the interaction 
is turned off after time t = Trj2 the amplitude of the den¬ 
sity oscillations would be enhanced. To illustrate this we 
assume that the coupling constant is a periodic function 
of time as in the inset of Fig. 3. We assume that the 
period is tt, i.e. the interaction gN = 0.3 is turned on at 
t = 0 then it is constant until t = nj2 where it changes 
sign and so forth. Obviously, in the homogeneous case 
the density of the impurity will not be altered. On the 
contrary, for the trapped case we should see that the 
amplitude of the density oscillations increases since the 
intrinsic time scale given by the harmonic trap coincides 
with the period for the coupling constant. To investigate 
this resonant behavior, we use the Green’s function from 
Eq. (37) with Eq. (38), which is easy to implement nu¬ 
merically for long time scales. The result of this investi¬ 
gation is shown in Fig. 3 and Fig. 4. In Fig. 3, we present 
|/(x,t)p at four different times, while Fig. 4 shows the 
evolution of the density from t = 0 to t = 197r/2. We 
clearly see that the amplitude of the oscillations at times 
equal to integer multiples of tt becomes much larger which 
means that higher natural orbits of the harmonic oscil¬ 
lator are populated. This demonstrates that even a very 
weak interparticle interaction appropriately modulated 
in time can drastically modify the initial wave function 
of the impurity. 



X 


FIG. 3. (Color online) Density profiles, \I(x, t)|^, as a function 
of the coordinate x taken at different times: t = 0 is shown 
by the solid (blue) line; t = 7r/2 - dashed (orange); t = 97r/2 
- dot-dashed (green); t = Stt - dotted (red). The inset shows 
the interparticle interaction which is a periodic function of 
time with period tt and peak amplitude gN = 0.3. t is in 
units of 1/D, X is in units of yjh/ (m/D), \I{x^t)\^ is in units 
of yj{rniQ)/h and g is in units of hTL^Jh/ (m/D) 



0.9 

0.6 

I 

0.3 

0 
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FIG. 4. (Color online) Density profiles, |/(x, t)|^, as a function 
of the coordinate x for different times t = /c7r/2 with k G 
{0,1, 2,..., 19}. t is in units of 1/D, x is in units of 
and |/(x,t)|^ is in units of {rniQ)/h. 


To quantify the evolution of the system we calculate 
the overlap between the states Ti(x) and /(x, t). That is 
we calculate 



Ti(x)/(x, t)dx 


2 


(41) 


From the definition it is obvious that Si{t) determines 
the decomposition coefficients of in the full ba¬ 
sis From another point of view So{t) compares 

the evolution of the interacting system with the non¬ 
interacting system and it is often referred to as the 
Loschmidt echo in the literature [57, 58]. Since the parity 
is conserved, Si{t) = 0 if i is an odd integer. We present 
Si{t) for i = 0,2,4 in Fig. 5. The Si for i > 6 are very 
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0 5 10 15 20 25 30 

t 

FIG. 5. (Color online) The overlap Si{t) from Eq. (41) as 
a function of time t. The (blue) circles correspond to i = 0, 
while the (yellow) triangles and (green) squares show the cases 
z = 2 and z = 4, respectively, t is in units of 1/Q. 



small for the time scales we consider. The figure clearly 
shows that the quench protocol from the inset in Fig. 3 
allows one to reach small values of the Loschmidt echo, 
So{t), and enhance 5'i>o even with weak interactions. 

For 2D and 3D a similar investigation for weakly- 
interacting systems can be provided. However, the dy¬ 
namics there can be understood from the ID case. In¬ 
deed the dynamics in the system should be driven with 
the Green’s function similar to Eq. (37) with additional 
complications arising from the degeneracy of the oscilla¬ 
tor levels. This complication can be lifted by noticing 
that angular momentum should be conserved. For ex¬ 
ample if the impurity is initially in the ground state, we 
expect the corresponding wave function to have zero an¬ 
gular momentum also for t > 0. From the same consider¬ 
ations as above, we conclude that the density oscillations 
for an interparticle interaction constant in time will hap¬ 
pen in the radial direction with period t = ir. Since this 
is very similar to the ID case, we do not go into more 
details here. 


V. SUMMARY AND OUTLOOK 

In this paper, we consider an ideal Bose gas with an 
impurity in harmonic trap. Assuming that the boson- 
impurity interaction is turned on at time t = 0, we inves¬ 
tigate real-time dynamics of the impurity. We present 
a field-theoretical approach to this problem using the 
ladder approximation. First, we show how to calculate 
the in-medium T-matrix in the space-time domain which 
form the basis for studying the dynamics of the system. 
As an explicit example, we address a weakly-interacting 
system in one spatial dimension. We build an infinites¬ 
imal time evolution operator to investigate the motion 
of the impurity. By repeated use of this operator, we 
show that the density of the impurity starts to oscillate 
with a period determined by the trapping frequency. We 
also show that the amplitude of the oscillations can be 
enhanced by periodically turning the interactions on and 
off. 


A number of interesting questions can be addressed 
directly using the ideas presented here. First, one can 
extend our approach to non-isotropic traps in order to 
address quasi-one (two)-dimensional setups which are of¬ 
ten used in experiments. The corresponding T-matrix 
can be easily found assuming different trapping frequen¬ 
cies in different directions. Second, while our formalism 
takes into account different masses for the bosons and 
the impurity, we have concentrated on equal masses for 
our numerical example. However, different masses would 
strongly affect the dynamics of the impurity. For in¬ 
stance, if the bosons are much lighter than the impurity 
then the impurity effectively feels a homogeneous distri¬ 
bution of the condensate. For the weakly-interacting case 
this would lead to absence of density oscillations. On the 
other hand, if the impurity is very light it should behave 
like a particle in a harmonic trap with the delta function 
peak generated by the bosons in the centre [59]. Other 
obvious extensions are the effects of stronger interactions 
and different initial states, i.e. when the impurity ini¬ 
tially has a density distribution different from the one 
that was considered. Naturally, the most interesting ex¬ 
tension is to go beyond the diagrams that were included 
here. For example three-body diagrams should lead to 
some manifestation of Efimov physics [60, 61]. Another 
interesting extension would be the investigation of the 
two-point correlation functions in the enviroment. These 
functions can be used to describe real time correlations 
between a boson and the impurity or two non-interacting 
impurities after the boson-impurity interaction quench. 

Most of the experiments with Bose gases are carried 
out with a weak boson-boson interaction. So a natural 
extension of the presented work would be to investigate 
the effect of this interaction on the time dynamics of 
impurity. Under the assumption that the profile of the 
condensate is not changed significantly by back action of 
the impurity one can apply the mean-field treatment to 
describe the ground state and low-lying excitations of the 
Bose gas [54] . The real time dynamics of the system then 
should lead to the depletion of bosons from the ground 
state neglected in this paper. 

Our study with an effective single-channel boson- 
impurity interaction of zero range applies to broad Fes- 
hbach resonances where a description of the interaction 
in terms of the scattering length is sufficient. However, 
we note that there are also some experiments with nar¬ 
row Feshbach resonances where the range effects can be 
significant [15], and therefore should be addressed in the 
future. 

Finally, we note that even one impurity can signifi¬ 
cantly modify the density of the Bose gas. For example, 
for a very heavy impurity placed in the middle of the trap, 
one can expect the depletion after time 1/Cl to be very 
large if the interaction is sufficiently strong. To describe 
this effect, one needs to carry out a more sophisticated 
calculation allowing many particles to be depleted from 
the condensate. Such a calculation would allow one to 
study the depletion of the condensate in real time. 
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Appendix A: Scattering length in 2D and 3D. 

Three-dimensional ease. For short-range potentials the 
zero-energy wave function at infinity has the form ^ (x) ^ 
1 — a/|x|, which is determined with only one parameter 
- the scattering length, a. One can establish [47] the 
following connection between the scattering length and 
potential: 


“^2^ / (^1) 

where 4^(x) is the wave function for the relative coordi¬ 
nate, 4^(x ^ oo) = 1 — If the potential is of short 
range we can approximate this wave function inside of 
the potential with constant, which allows us to rewrite 
this definition as 

where |xe| is the range of the interaction. As shown 

I 


below, for the Gaussian potential considered in the text 
|xe| = c/A with c = ^/n/2 and 


a = 


9h 

2TTh? 



(A3) 


Notice that for any non-zero value of scattering length, 
the limit of A ^ oo corresponds to ^ ^ 0 . 

Two-dimensional ease. In two spatial dimensions the 
standard way to define the scattering length, 020 > 0 , 
for short-range potentials is also through the zero-energy 
wave function of relative motion, which at infinity should 
be of the form ln(|x|/a 2 D) [62]. This form can be ob¬ 
tained directly by solving two-body problem, see for ex¬ 
ample [63], from which one can deduce 


\n{a2D/l) 


-1 + /q dssln{s/l)V{s)^{s) 
^ dssV{s}^{s) 


(A4) 


with 4^(0) = 1 and I any parameter with dimension of 
length. Again assuming that the interaction is of short 
range, we put = 1 , which for the Gaussian potential 
leads to the following dependence of scattering length on 
g and cutoff parameter A 


^2 

ln(a 2 DA) = (A5) 

Z jig 

where 7 = 0.57721... is the Fuler-Mascheroni constant. 
Similar to the 3D case A ^ oc corresponds to ^ ^ 0. 


Appendix B: Regularization procedure for the T-matrix 


In this Appendix, we show how Eq. (20) is obtained. We also present the corresponding regularization procedure in 
2D and 3D. Noting that /(t) is independent of the spatial coordinates, we show the derivation for the homogeneous 
case, i.e. (2 = 0, which is more transparent. We first perform a transformation to the momentum-time domain, since 
the Green’s functions conserve momentum in the homogeneous case, i.e. 

G{p,t,p',t') = {2'K)’^G{p,t,t')S{p - p'). (Bl) 


Moreover, we define the T-matrix in the momentum-time domain as 


^{Pi,Pi,ti,P2,P2,t2) = J ^{xi,x[,ti,X2,X2,t2)e (B2) 

where p denotes a wave vector. With this definition, Eq. (9) assumes the following form 

^iPl,Pl,h,P2,P2,t2) =V{P2 -Pl)5{ti -t2)5{pi +p[ -P2 -p' 2 ) 

J ^(Pi ~ r)G^j{p\ +pi — K,,ti,k',t)G%{K,,ti, k!, t)T{k', k!, t,p2,p'2,t2)dk'dndn!At. 


h{27r) 


(B3) 
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We rewrite this equation using the Green’s function in free space from Eq. (Bl) 

^iPl,Pl,h,P2,P2,t2) =V{P2 -Pl)5{ti -t2)Sipi+p'i -P2 -p' 2 ) 

y y(pi - k)G°{p\ +Pi- +pi- K,K,t,P2,P2,t2)dKdt. 


h{27r) 


D 


(B4) 


Furthermore, we introduce a new set of variables P = pi p'l, K = p 2 -\-P 2 ^ 2 p = pi —pi^2k = P 2 — P 2 ' Since the total 
momentum is conserved, we define new function T from r(pi,p'^, ^ 2 ) = ^{P — ^)r(P;p, /c, ^ 1 ,^ 2 ) which leads 

to 


f(P;p,/c,ti,t2) =V{k - p)S{ti - t2) 


+ 


h{27r) 


D 


j E(P/2 — p —/i:)G 5(P —ti, t)G^(/^, ti, t)f(P; P/2 —/^,/c, t, t2)d/^dt. (B5) 


To illustrate appearance of a pathology, we take V{k) = g where ^ is a constant satisfying 0 < ^ < 00 . For this 
potential, the T-matrix reads 


r{P;p,k,ti,t2) = gS{ti - t2) + 9 J ~ i^,ti,t)G%{K,ti,t)r{P; P /2 - K,k,t,t2)dKdt. 


(B6) 


Notice that F is independent of its second and third arguments. The integral over k can be taken easily, since the 
Green’s functions are just exponents, so 


t{P-,ti,t2) = gS{ti-t2) - g , b [ dtt{P;t,t2) 

[Zw) n Jti 


27r/i 


ih{t — ti) 


D/2 




(B7) 


The dependence on P can be eliminated easily, assuming the following form of the T-matrix r(P;ti,t 2 ) = 
g * 2 (mj+ms) ti )— ti), which leads to the equation for / 


f(T) = [ Orf(T (B8) 

This equation is presented and discussed in the main text. However, the equation is pathological. For example, the 
reasonable assumption that / is a non-zero function that is analytic almost everywhere cannot be satisfied as the 
integral has divergence at small times. 

Now we use the Gaussian potential defined in Eq. (18) and take the limit A ^ cxo afterwards. Using the results 
from above, we assume that T does not depend on its second and third arguments. This yields 


t{P-, 0 ,T) =g(^S{T) + dtt{P-,T-t) Ie-^G°j{P - K, 0 ,t)G%{K, 0 ,t)dK 


T / \ 


h^P^ f2 
4m ? 


.hP‘^ 


1 + 2 m/ 

4A2 ' 2/x ^ 


(B9) 


We can neglect value of 1/4 in the denominator of the argument in the exponent, since it contributes only for very 
small values of t and exponent is a bounded function. The factor of 1/4 in the denominator of the fraction cannot 
be neglected as (1/t)^/^ can be an arbitrarily large number for small t. As before we assume that r(P;ti,t 2 ) = 
e ^ 2 (mj+ms) ^^2 ti) Icads to 






(BIO) 


where we extended the integration limits to infinity since the function / should vanish for negative values of the 
argument. Notice that assumption that / is analytic almost everywhere and non-zero can be satisfied now since the 
divergence of integral for A ^ oc is compensated by the coupling constant which goes to zero. It is simpler to 
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take the limit of zero-range interaction after the expression is Fourier transformed to the frequency domain, where we 
obtain 


The integral can be carried out and has the same kind of divergence as l/^f, so the limit A ^ oo yields a finite 
expression that depends only on the scattering length. For instance in 3D, we have 

(1 - aA/c)/(w) = f{u!){-SiiJ,A + 4\/27rw/i^/^)^. (B12) 

This equation can be used to determine c which eliminates the divergence proportional to A so c = y/TTl2. After this 
divergence is eliminated we are left with the equation for /(cj), presented in the main text. In 2D, we have instead 

-^(ln(a2^A) + 7/2 )/(u;) = 1 + (in (^) + 1 - in) , (B13) 

which yields the result shown above. 


Appendix C: r(0;0,t) in the one-dimensional case 


To obtain f{T) in the main text, we applied the inverse Fourier transform. In this appendix, we show another way to 
derive /(T). We focus on ID, since there regularization of the integrals is not needed, by considering f{t) = r(0; 0, t), 
which is written as 


f (0; 0, T) =g (^6{T) + dtf (0; T - 0, 0, 


(Cl) 


By definition this integral is understood as the following sum 

n pT pT pT 


/ f in \ ^ pT n p \ 

f(0;0,T) =5((5(T) + y] (—j J dti J dt2... J dtn-i'[[jdKjG°{Kj,tj-i,tj)GB{Kj,tj-i,tj)], (C2) 

n=l ^ tn 2 (/ = 1 

where to = 0,tn = T; for n = 1 it is implied that the integral over time is absent. The integrals over Hj are easily 
taken, such that 


f(0;0,T) =flp(r) + 




1 


i^j ^i-i) 


(C3) 


From dimensional analysis the multidimensional integral should be equal to anT'^^‘^ To find the coefficients, an, 
we establish the recursive relation 


rpnj 2 — 1 _ 

n.-L — 


pT pT pT n I I pT 1 « 

dti dt2.: ----=/ -^an-i{T-ti) 

Jo Jtl Jtr,-2 Jo Vh 


(C4) 

and obtain Since ai = I we obtain a 2 k = (k^-i)\ ^ <^ 2 /c+i = (/g-i)! large 

values of k we have B{l/2^k) ^/n/k. This allows us to conclude that the coefficients for large k decay roughly 

as 7r’^/^/(n/2 — I)!, which implies that the sum in Eq. (C3) is convergent for any g and T > 0. Finally this sum is 
written as 


\ n=l \ / \ / 


n/2 


(C5) 


We notice that this expression for ^ > 0 converges to Eq. (27). Eor 2D and 3D similar ideas can be used, however, a 
proper regularization procedure is needed which makes the calculations more complicated. 
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Appendix D: Molecular branch in 3D. 

In 3D a > 0 implies existence of a bound molecule, which then smootly evolves to a Bose polaron [46]. This smooth 
evolution should be contrasted with the Fermi polaron case [27]. In the frequency domain the molecule is seen through 
the pole close to the integration contour, which leads to existence of a double solution of the standard equation for 
the energy of a quasiparticle 


uj = nRe/(co’). (Dl) 

It is interesting to note that the double solutions also exist in ID for ^ < 0, where Auj = 

— [g‘^jj. ± — Sn — Agn). Note that in this section we put h = 1 for simplicity. We also note that, in 

ID for ^ ^ 0 Eq. (^Dl) not always can be used, at least within the T-matrix approximation, to define a quasiparticle 
due to strong correlations happening in ID. In this section we would like to investigate appearence of this molecular 
solution in 3D for small postitive a in the time domain. For this we write the equation for / for /i = I using the 
derivations from Appendix B 

/ 1 \ ^/2 /■* / 1 

= I (—^) (D2) 

For very large times this equation reads 

/ 1 \ 3/2 nt^oo / 1 \ 3/2 

with one obvious solution f{t oo) = that describes propagation of a molecule. So we look for / in the form 
f{t) = K{t)e^ ^ where K(t) solves the following equation 

/ -j \ 3/2 pt / -I \ 3/2 

K{t)=g5{t)-igi^—j (D4) 

Let us assume that K is constant. With this assumption we obtain 

K = 2„.(r)-^ 

where erf(x) is the error function. The equation is satisfied for small a if A = —inf a. This means that in the limit of 
small scattering length propagation of the virtual dimer is determined by the following function f{t) = ^ 

2na5{t). This result is very obvious, since it also should be reached if the interparticle interaction is completely 
repulsive, such that a molecule cannot be formed. 
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